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fSJ ' Abstract 

(N ■ 

We study chiral symmetry structure at finite density and temperature in the presence 
of external magnetic field and gravity, a situation relevant in the early Universe and in the 
' core of compact stars. We then investigate the dynamical evolution of phase transition in 

. the expanding early Universe and possible formation of quark nuggets and their survival. 

(N 
O 

CN ■ 1 Introduction 
O 

Q^' Spontaneous symmetry breaking is one of the most important concepts of all unified gauge 
theories. The idea that underlying symmetries of nature are larger than that of the vacuum 
D ! play a crucial role in the unification of forces. Of particular interest is the expectastion that 
' at high temperatures, symmetries that are spontaneously broken today are restored and that 
during evolution, the Universe passed through a series of phase transitions from a higher 
^ . symmetric phase to a lovi^er symmetric phase associated with the spontaneous breakdown of 
^ \ gauge or perhaps global symmetry. In particular at t = 10~^'^s when the temperature was 
~ 2Q{)GeV , the Universe passed through the EWS breaking phase transition and later at 
t = 10~^s at T = 200MeU, there must have been a QCD phase transition from QGP to 
confined hadronic matter and also to chiral symmetry breaking phase transition. Since the 
vacuum structure of SSB theories is very rich, topologically stable configurations of gauge 
and Higgs fields in the form of domain walls, cosmic strings and monopoles on the one hand 
and non-topological solitons like Q balls, quark nuggets, solitone stars on the other may exist 
and may have observable signature. It is thus very instructive to investigate how the phase 
transition takes place in QFT in the environment of the early Universe and in the core of 
neutron stars where temperature, density, external electromagnetic field and external gravity 
may all play important role. 

In GUTS Higgs play a most important role. They are fundamental scalars which give masses 
to fermions and gauge bosons through their VEVs. In theories like Technicolour models, 
Higgs are composite of some fundamental fermion fields. QCD is an example of QFT which 
is invarient under chiral transformations at the Lagrangian level in the absence of quark mass 
matrix, the dynamics of QCD are expected to be such that chiral symmetry is dynamically 
broken with the vacuum state acquiring a non-zero quark-antiquark condensate < qq >, and 
the Goldstone theorem then requires the existence of approximately massless pseudoscalar 
mesons. To study chiral phase transition in QCD we need a non-perturbative treatment and 
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a particularly attractive frame work to study is the Nambu-Jona-Lasino model. The linear 
Sigma model is another such model which has the advantage of being a renormalizable model 
in contrast to the NJL model which in 4-D is known to be non-renormalizable. An elegent 
and efficient way to study symmetry properties of vacuum at finite temperature and density 
in external envirornment is through the 'Effective Potential' approach discussed extensively in 
the litrature. We compute here, in the one loop approximation, the effective potential in the 
presence of external magnetic field and gravity at finite temperature and density in the frame 
work of linear sigma and NJL model respectively. We then study the dynamic evolution of 
phase transition through bubble nucleation of hadronic phase. 



2 Chiral symmetry in external magnetic field 

It has also been suggested that systems with spontaneously broken symmetries may make 
a transition from broken symmetric to restored symmetric phase in the presence of external 
fields. Large magnetic fields with strength upto 10^^ gauss have been conceived to exist at 
the time of supernova collapse inside neutron stars and in other astrophysical compact objects 
and in the early Universe. Effect of such a strong magnetic field on chiral phase transition is 
thus of great interest for baryon free quark matter in the early universe and for high density 
baryon matter in the core of neutron stars. To study chiral phase transition in QCD we need 
a nonperturbative treatment. Lattice techniques and the Schwinger-Dyson equations provide 
specially powerful methods to study the chiral structure of QCD. A particularly attractive 
frame work to study such systems is the linear sigma model originally proposed as a model for 
strong nuclear interactions. We will consider this as an effective model for low energy phase 
of QCD and will examine the chiral symmetry properties at finite density and in the presence 
of external magnetic field. To fix ideas we consider a two fiavor SU (2) x SU (2) chiral quark 
model given by the lagrangian. 

£ = iV^T^a^V- - 9^{(^ + i75T.7r)V' + ^{df^a)^ + \{d^^f - U{a,^) (1) 

where ip is the quark field a and vf are the set of four scalar fields and g is the quark meson 
coupling constant. The potential V{a, vr) is given by 

C/(a,vr) = -i/x2(a2 + vf2) + iA(a2 + vfY (2) 

For /J? > chiral symmetry is spontaneously broken. The a field can be used to represent 
the quark condensate, the order parameter for chiral phase transition and the pions are the 
Goldstone bosons. At the tree level the sigma , pion and the quark masses are given by 

ml = ZXali - i?;ml = Xali - //^ = ga^i (3) 

where a^i = ^ = f^. We compute here, in the one loop approximation, the effective potential 
in the presence of external magnetic field which is defined through an effective action T(a, B) 
which is the generating functional of the one particle irreducible graphs. The effective potential 
is then given by 

Veff{a,B) = Vo{a) + Vi{a,B) (4) 

where Vi{a, B) is obtained from the propagator function G{a, B) by the usual relation Vi{a, B) = 
— ^Tr log G{a, B). Alternatively one can compute the shift in the vacuum energy density due 
to zero-point oscillations of the fields considered as an ensemble of harmonic oscillators |^]. 
We thus require energy eigenvalues (excitations) of particles in the magnetic field, which can be 



easily obtained, and in the absence of anomalous magnetic moment for uniform static magnetic 
field in the z-direction for a particle of mass M, charge q and spin J, are given by 



E{k,,n, J,) = {kl + M2 + (2n + 1 - sign{q) j,) \q\ sf (5) 

where n represents the landau level. In the presence of magnetic field, all we need to do is to 
replace the phase space integral / by |^ J2'^=o (2np ^^'^ energy by expression (|5|) for 
charged particles. The one loop effectiv potential can now be calculated by standard techniques 
by dimensional regularization. The effective potential has poles which can be absorbed in 
counter terms. The finite part depends on the exact renormalization conditions that are 
imposed. In what follows we would use the MS renormalization scheme. To proceed further 
we first consider the case < 1 , keeping leading terms and adding the contributions of all 
the charged particles the total V(,ff{a,B) for the sigma model at the one loop level is thus 
given by [§ 

Veff{a,B) = -\^^a^ + \a' 
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For > 1 ,the last term in eqn.(^) is replaced by ^'^^g^ [l — logM"^) summed over fiavors. In 
figure 1, we plot Veff{a, B) as a function of a for different values of magnetic field and compare 
it with the case of zero magnetic field by. As input parameters we choose the constituent 
quark mass mf = 500 MeV, sigma mass mo-=1.2 GeV and /7r=93 MeV. We find that in the 
presence of intense magnetic fields the chiral symmetry breaking is enhanced. For magnetic 
field large compared to m^, we observe that though the fermionic contribution is towards 
symmetry restoration, it is not enough to offset the contribution of charged goldstone pions. 
In order to study chiral symmetry restoration in the case of neutron stars as a function of 
chemical potential /U associated with finite baryon number density we employ the imaginary 
time formalism by summing over Matsubara frequencies. This amounts to adding the fermionic 
free energy to the one loop effective potential and is given by 

which in the presence of static uniform magnetic field becomes 



where 7 is the degeneracy factor and is equal to 2N(. for each quark fiavor. We consider cold 
dense isospin symmetric quark matter for which the integrals can be performed analytically. 



The baryon number density corresponding to the chemical potential /i is given by the usual 
thermodynamical relations. 
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for zero and finite magnetic field, respectively. Here Timax — 

Int r2\q \ B • To study chiral 
symmetry behavior at finite density in the presence of uniform magnetic field, we minimize 
effective potential with respect to the order parameter a for fixed values of chemical potential 
and magnetic field ( which then fixes the baryon density ). The results are shown in figure 
2 where we have plotted the order parameter cr as a function of density at T=0 for different 
values of magnetic field. The solution indicates a first order phase transition. The actual 
transition takes place at the point where the two minima of the effective potential at cr=0 and 
a=a(ii,B) non zero become degenerate. The lower values of a (shown by dotted curves) are 
unphysical in the sense that they do not correspond to the lowest state of energy. We find 
that magnetic field continues to enhance chiral symmetry breaking at low densities as expected 
but as the magnetic field is raised the chiral symmetry is restored at a much lower density 
compared to the free field finite density case. This can be clearly seen from figure 3 where we 
have plotted the phase diagram in terms of baryon density and magnetic field. 



3 Chiral symmetry in NJL model in curved space time 

The NJL model four fermion theory in fiat space-time in arbitrary dimensions has been studied 
1^] using the 1 /N expansion method. It is shown that chiral symmetry is restored in the theory 
under consideration for sufficiently high temperature or chemical potential. It is found that for 
space-time dimensions 2 < D < 3 both a first order and second order phase transition occur 
depending on the value of the four fermion coupling while for 3 < D < 4 only the second order 
phase transition exists. Here we study the chiral phase structure of the four fermion theory 
in curved space-time at finite temperature and density. The Nambu Jona Lasinio model in 
curved space time is defined by the action 

S = y'ci^xy(^[i^7^(x)VM^ + ^( (V5V')' + (V5i75rV')')] (H) 

where g is the determinant of the space time metric ,7^(2;) the Dirac matrix in curved space- 
time, VajV' the covariant derivative of the fermion field ip and N is the number of colours, 
we take the number of flavors to be two. We work in the scheme of the 1/N expansion and 
perform our calculations in the leading order of the expansion. For practical purposes it is 
more convenient to introduce the auxiliary fields a and vf and consider the equivalent action. 

S = J d^x^i^) [i i^rix) V/. V' - ^(<t' + ^') -i^ia + i 75 r 7?)^] (12) 

Replacing a and vr by the solutions of the Euler-Lagrangian equations arising from (12) we 
reproduce the action (11). This is because the fields a and 7? are not independent degrees of 
freedom in (12) and the Euler-Lagrangion equations for a and tt are infact constraint equations 
which fix a and vf given tp and ■0. If a non vanishing vacuum expectation value is assigned to 



the auxiliary field a, then there appears a mass term for the fermion field ip and the discrete 
chiral symmetry is eventually broken. The effective potential (with N factored out) in the 
leading order of the 1 /N expansion is then given by 



V{a, vr) = — (f^ + 71"^) + iTrln S {x, x; s) 
2 A 



(13) 

=(T+ 475 if 

In equation (13) the variables a and vf are regarded as constantand the Green function S is 
the solution of the equation 

(i7^(x) Vm -s)S{x, y; s) = —l==d^{x - y) (14) 

v(-f(2;)) 

Thus the effective potential is described by the two point Green's function S(x,x;s) of the mas- 
sive free fermion in curved space-time. Using the Green function obtained in the approximation 
of keeping only linear terms in the curvature, the effective potential which in D-dimensions 
read as follows: 

'2 '•'^ I- d^k 1 
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where J"'' = ^[7'*, 7^] and latin indices are vierbein indices. The effective potential is divergent 
in two and four dimensions and is finite in three dimensions in the leading 1/N expansion. The 
four fermion theory is renormalizable in 2-D flat space. Using the renormalization condition 

^ where M is the renormalization scale. The renormalised effective potential 
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in D dimensions is given by 
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xr(i-y)^ (16) 

We now obtain the four dimensional limit of the NJL model in the MS renormalization scheme 
given by 

V(a, 0) 1,0", 9 1, x/<7,9 

^ ( /n(^)2 - 1 -/n 4^ + 7) (^)2 (17) 



Alternatively one could regularize the divergent part by cutting off the momentum integral at 
finite cutoff A [0. This gives 



' --'lnil + ^) + J^] (18) 
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The two expression can be shown to be equivalent after carrying out the renormahzation of 
the couphng constant. The ground state of the theory is determined by the minimum of the 
effective potential (17) namely, by solving the gap equation 



=0 (19) 
da ''^='" ^ ' 



For Xr > Ac = i-j_3;„4 7r-3 7 i minimum of the effective potential is located at the non- 
vanishing cT, the chiral symmetry is broken down dynamically and a dynamical fermion mass 
is generated. At the critical point the effective potential has two degenerate local minima 
obtained by putting 

V{a,0) = V{m,0) = (20) 

The solution of the gap equation and the value of the critical curvature can be obtained numer- 
ically. First, we fix the coupling constant A,- greater than the critical value Ac corresponding to 
the broken symmetric phase. To study the phase structure in curved space-time we evaluate 
the effective potential (17) by varying the space-time curvature. We see from Fig. 4 that the 
chiral symmetry is restored as R is increased for a fixed A. The phase transition induced by 
curvature is of first order. The effective potential V^{a,0) at finite temperature, density and 
curvature can be obtained by standard techniques as in Section 2 and is given by 

V^{a,0) = (l + 3/n4^-37) ^2 

-8^^'^M^-2-^^^^ + ^)^ 
R 

'' T79 — ^ —InA-K + 7)0"^ 
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At high densities and low temperatures or at high temperatures and low densities, the integrals 
in (21) can be carried out analytically by the standard techniques. In Fig. 5 we have shown 
the phase boundary in temperature, chemical potential plane for different values of curvature. 
We notice that with the increase in R, phase transition takes place at lower temperature and 
density. In Fig. 6 we have shown the phase boundary curve in the temperature, curvature and 
chemical potential plane. 



4 Dynamical evolution of quark-hadron phase transition through 
bubble nucleation 

It is well known that a phase transition from quark gluon plasma to confined hadronic matter 
must have occurred at some point in the evolution of the early Universe, typically at around 
10 — 50/is after the Big Bang. This leads to an exciting possibility of the formation of quark 
nuggets through the cosmic separation of phases . As the temperature of the Universe falls 
below the critical temperature Tc of the phase transition, the quark gluon plasma super cools 



and the transition proceeds through the bubble nucleation of the hadron phase. The typical 
distance between the nucleated bubbles introduces a new distance scale to the Universe which 
depends critically on the super cooling that takes place. As the hadronic bubbles expand, they 
heat the surrounding plasma, shutting off further nucleation and the two phases coexist in 
thermal equilibrium. The hadron phase expands driving the deconfined quark phase into small 
regions of space and it may happen that the process stops after the quarks reach sufficiently 
high density to provide enough pressure to balance the surface tension and the pressure of the 
hadron phase. The quark matter trapped in these regions constitute the quark nuggets. The 
number of particles trapped in the quark nugget, its size and formation time are dependent 
sensitively on the degree of super cooling. The duration of the phase transition also depends 
on the expansion of the Universe and on other parameters like the bag pressure B and the 
surface tension a. 

The quark nuggets formed in the small super cooling scenario are in a hot environment 
around the critical temperature Tc and are susceptible to evaporation from the surface |^] and 
to boiling through subsequent hadronic bubble nucleation inside the nuggets [0]. However in 
the large super cooling scenario we have the interesting possibility of these nuggets forming at 
a much lower temperature than Tc due to the long duration of the transition and consequent 
expansion of the Universe. Alcock and Farhi have shown that the quark nuggets with 
baryon numbers < 10^^ and mass < IO^^Mq are unlikely to survive evaporation of hadrons 
from the surface. Boiling was shown to be even more efficient mechanism of nugget destruction 
|To| . These results were somewhat modified by Madsen et.al. jll|] by taking into account the 
flavor equilibrium near the nugget surface for the case of evaporation and by considering the 
effect of interactions in the hadronic gas for the case of boiling. In the large super cooling 
scenario the time of formation of these nuggets can be quite late when the number of baryons 
in the horizon (of size ~ 2t) is large and temperature much lower. These nuggets can easily 
survive till the present epoch. 

There have been recent observations by gravitational micro lensing [|l2| of dark objects 
in our galactic halo having masses of about 0.01 — 1 solar mass. If these objects have to be 
identified with quark nuggets, they could only have been formed at a time later than the time 
(~ 50 — lOO^s) when the Universe cooled through Tc. S When the early Universe as a quark- 
gluon thermodynamic system cools through the critical temperature Tc, energetically, the new 
phase remains unfavorable as there is free energy associated with surface of separation between 
the phases. Small volumes of new phase are thus unfavorable and all nucleated bubbles with 
radii less than critical radius collapse and die out. But those with radii greater than the 
critical radius expand until they coalesce with each other. So super cooling occurs before the 
new phase actually appears and is then followed by reheating due to release of latent heat. 
The bubble nucleation rate |jl^ at temperature T is given by 

Wr 

I = he-— (22) 

where Iq is the prefactor having dimension of T^. The prefactor used traditionally in early 
Universe studies is given by Iq = (^^)^^^- Csernai and Kapusta 10 have recently 
computed this prefactor in a coarse-grained effective field theory approximation to QCD and 
give lo = ^(^)^ g^(Aw)^ where ??g = lAAT^ is the shear viscosity in the plasma phase, 
is a correlation length of order 0.7 fm in the plasma phase and AtD is the difference in the 
enthalpy densities of the two phases. In this letter we use both these prefactors for comparison. 
The critical bubble radius Rc and the critical free energy Wc are obtained by maximizing the 
thermodynamic work expended to create a bubble and are given by Rc = p^^^j^^^p and 

Wc = 47rcTi?2(T)/3 = where AP = Ph{T) - Pg{T) is the pressure difference in hadron 

and quark phase. 



For simplicity we describe the quark matter by a plasma of massless u,d quarks and 
massless gluons without interaction. The long range non-perturbative effects are parameterized 
by the bag constant B. The pressure in the QGP phase is given by Pq{T) = ^gq^T^ — B 
where Qq ~ 51.25 is the effective number of degrees of freedom. In the hadronic phase the 
pressure is given by Ph{T) = ^gh^T"^ where gh ~ 17.25, taking the three pions as massless. 

The fraction of the Universe h{t) which has been converted to hadronic phase at the time 
t is given by the kinetic equation 



h{t)= fl{T{t'))[l-h{t')]V{t',t) 



^^^'^^'dt' (23) 
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where V(t', t) is the volume of a bubble at time t which was nucleated at an earlier time t' and 
R{t) is the scale factor. This takes bubble growth into account and can be given simply as 



Rcim) + £ ^MT{t"))dt" 



(24) 



where v{T) is the speed of the growing bubble wall and can be taken to be v{T) = — t^I] 2 
where Vq = 3c. This has the correct behavior in that closer T is to Tc slower do the bubbles 
grow. When T = we have v{T) = ^ the speed of sound of a massless gas. For T < 

which occurs when there is large super cooling, we use the value v{T) = 

The other equation we need is the dynamical equation which couples time evolution of 
temperature to the hadron fraction h{t). We use the two Einstein's equations as applied to 
the early Universe neglecting curvature. 



R /SttG 1 

^ = -^^ (26) 
R 3wdt ^ ^ 

where w = p + P is the enthalpy density of the Universe at time t. The energy density in the 

mixed phase is given by p{T) = h{t)ph{T) + [1 — h{t)]pq{T), where ph and pq are the energy 

densities in the two phases at temperature T and similarly for enthalpy. We numerically 

integrate the coupled dynamical equations (23), (25) and(26) to study the evolution of the 

phase transition starting above Tc at some temperature T corresponding to time t obtained by 

integrating the Einstein's equation (25) and (26) in the quark phase. The number density of 

nucleated sites at time t, is given by 

Nit) = J^'l{t')[l-hit')](^^ydt' (27) 

Therefore the typical separation between nucleation sites is Z„ = N{t)~^^^. This distance scale 
will eventually determine the number of quarks in a nugget. This scale can be up to 10^^ Km 
depending on the parameters B and a which correspond to a distance of ~ lAMpc today. 
The nuggets could be points in space around which, later in time, the matter in the Universe 
may have gravitationally clustered to give the observed large scale structures in the Universe. 
The observable separation of galaxies in the Universe can be a remnant of this transition with 
the centers of the galaxies being the quark nuggets. Of course the collision of bubbles and 
their random nucleation and interaction will also lead to clustering of the nuggets, which can 
qualitatively explain the clustering of galaxies. To get an idea about the super cooling before 



nucleation begins, we can plot the nucleation time as a function of temperature, defined by 
Enucleation = The quark number density is given by Ug = ^C{3){^)T^ where ^ is 

the quark to photon ratio estimated from the abundance of luminous matter in the Universe 
to be roughly equal to 3 x 10^^'^. The quark nuggets have Nq quarks at time t given by 
the number of quarks in a volume i-e. Nq = 7^^- The nucleation sites are actually 

randomly distributed, but we expect a distribution of quark numbers around Nq. The average 
temperature at which nuggets are formed when bubbles coalesce is obtained by finding the 
average time at which the expanding bubble surfaces meet. Assuming a cubic lattice, we have 
done this numerically to get the corresponding time tf and temperature Tf for different values 
of B and a. When the fraction of the space occupied by the bubbles is around 50 percent, we 
expect the bubbles to meet in an ideal picture,i.e. if all bubbles are essentially nucleated at one 
instant which is the maximum nucleation time and they all have the same radius. However 
we have a distribution of expanding bubble sizes because of the different points of time at 
which they were nucleated. Therefore the estimate of the time of nugget formation by treating 
all bubbles to be of the same size is an underestimate. We find that hadron fraction h{t) is 
only around .12 when bubbles meet by this criteria. However we do not expect this to change 
qualitatively the broad picture of the transition and the nugget formation apart from reducing 
the formation time. 



Fig. 7 [15 1 shows the temperature as a function of time. It is clear from this diagram that 
reheating takes place as nucleation starts with the release of latent heat. As a increases, the 
super cooling is larger and reheating is slower. The transition takes much longer to complete 
with more chance of nugget formation. This allows the nuggets to be formed at a much 
lower temperature when bubble walls meet. For low supercooling there is rapid reheating, 
temperature reaches Tc and the phase transition is completed very swiftly with no chance of 
any nugget formation. Fig. 8 shows the fraction h{t) of the Universe in hadron phase as a 
function of time. For small values of a the transition completes quickly as h{t) goes to 1. 
But for larger a it takes a larger time for h{t) to become 1. We also notice that in the large 
super cooling scenario the Kapusta prefactor becomes much bigger than the standard one by 
many orders of magnitude. This makes the nucleation rate as well as the reheating faster. In 
the case of low super cooling the two prefactors give almost identical results. The number of 
quarks in the horizon NqH at time t is Ngfj ~ ng(|7rt^) ~ i'fr) '^^^^ T^Ant^ and we find that for 
all interesting cases Nq < NqH and this number is very sensitively dependent on the surface 
tension. Physically it is possible to have Ngjj > Nq > 10^^ for some values of the parameters 
B and a. In table I below we list some physical quantities for some representative values of B 
and a. 





a 


Tc 


tf 


Tf 




NgH 


^n 


MeV 


MeVfm-'^ 


MeV 




MeV 






m 


235 


50 


169 


12.1 


169 


2.6 X 10^^ 


7 X 10^^ 


8 X 10"^ 


145 


57.1 


104.4 


34.7 


99.9 


9.4 X 10^^ 


3.4 X 10^2 


4.6 


125 


57.1 


90 


56.1 


78.8 


1.1 X 10^9 


7.3 X 10^3 


63 


125 


77 


90 


1511.8 


17.9 


6.3 X lO''^ 


1.6 X 10^^ 


4.9 X 10^ 


100 


39.5 


71.9 


2595 


13.9 


1.2 X 10^° 


3.8 X 10^6 


8.4 X 10^ 


113 


57.1 


81.3 


5138 


12.3 


9 X 10^9 


2 X 10^^ 


1.7 X 10^ 



5 Conclusions 



In conclusion we have examined the chiral symmetry behavior of the Linear Sigma model in the 
presence of static, uniform magnetic field at the one loop level at zero density and at densities 
relevant in the core of neutron stars. We find that the contribution of scalar and fermion loops 
leads to an increase in chiral symmetry breaking. At high densities too, this effect persists 
and for magnetic fields of strength upto 10^^ Gauss, there is enhancement in chiral symmetry 
breaking resulting in the restoration of symmetry at densities higher than if no magnetic field 
were present. However, in the case of high magnetic field B > 10^^ Gauss the chiral symmetry is 
restored at lower densities. Thus in the core of neutron stars, if the nuclear matter undergoes a 
transition to deconfined quark matter, the presence of magnetic field would imply the existence 
of massive quark matter due to enhancement in chiral symmetry breaking. This would affect 
the equation of state and will have astrophysical implications. We then investigated the phase 
structure of NJL model in curved space time in the linear curvature approximation and found 
that in the presence of external gravity with positive R, the chiral symmetry restoration is first 
order even as in flat space-time, the transition is second order with temperature and becomes 
first order in the presence of density. We find that with the increase in R phase transition takes 
place at lower temperature and density. Regarding evolution of quark-hadron phase transition 
detailed dynamics of the quark-hadron transition in the early Universe show that the evolution 
of the Universe does not necessarily follow the small super cooling scenario and certain choices 
of B and a can have a bearing on the present state of the Universe. As nuggets with Nq > 10^^ 
are expected to survive the transition, they will contribute to the density of the Universe. We 
have explored in detail the possibility of nugget formation and also estimated their average 
separation, time of formation, quark content and survivability. Clearly, the analysis can be 
improved by taking interactions into account in both the phases and also bubble interactions 
may be incorporated in the calculations, the qualitative results however arc not expected to 
change. Thus if the nuggets studied above arc indeed formed in a much cooler environment, 
they could contribute significantly to the missing mass in the Universe and be candidates for 
dark matter. 
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viding the opportunity to present this work and for making the conference stimulating and 
enjoyable. 



Figure Captions 



Figure 1. Effective Potential in units of (lOOMey)^ as a function of (j{MeV) for different 
values of the magnetic field. The curves a, b, c, d and e are for B=0, 10^^,10^^,10^^ and 
3 X 10^^ Gauss respectively. 

Figure 2. Chiral condensate a{MeV) as a function of baryon density in for magnetic 
field B=0,10^6,10^^andlO^'^ Gauss respectively. 

Figure 3. Phase diagram as a plot of magnetic field versus baryon density in units of fm~^. 

Figure 4. Behaviour of Effective Potential as a function of jj for different values of 
the curvature for A > Acr(A = 10). The curves a, b, c, d, e and f are for R = 16, 13, 12, 8, 
and -2 respectively. 

Figure 5. Temperature vs chemical pot. for different values of curvature. The curves a, b, 
c, d, e and f are for R/M^ = 0, 2, 4, 6, 8 and 10 respectively. 

Figure 6. Phase diagram in the temperature, curvature and chemical potential plane. 

Figure 7. Temperature as a function of time. Solid and dashed curves are for = 
lOOMey and a = 39.5MeV fm~'^ with the standard and the Kapusta prcfactors respectively. 
Long dashed and dotted curves are for 5^/4 = llSMeV and a = 57.1MeV fm.-'^ with the 
standard and the Kapusta prefactors respectively. 

Figure 8. The hadron fraction as a function of time. Curves as in fig. 7. 

Table Caption 

Tabic I. Some relevant physical quantities for some representative values of B and a using 
the standard prefactor. 
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